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We present two-dimensional resistive magnetohydro dynamic simulations of line-tied 
asymmetric magnetic reconnection in the context of solar flare and coronal mass ejection 
current sheets. The reconnection process is made asymmetric along the inflow direc- 
tion by allowing the initial upstream magnetic field strengths and densities to differ, 
and along the outflow direction by placing the initial perturbation near a conducting 
wall boundary that represents the photosphere. When the upstream magnetic fields are 



^ il asymmetric, the post-flare loop structure is distorted into a characteristic skewed candle 

Q ' flame shape. The simulations can thus be used to provide constraints on the reconnec- 

^ ■ tion asymmetry in post-flare loops. More hard X-ray emission is expected to occur at 

j^ . the footpoint on the weak magnetic field side because energetic particles are more likely 

to escape the magnetic mirror there than at the strong magnetic field footpoint. The 
CN ! footpoint on the weak magnetic field side is predicted to move more quickly because 

JI^ I of the requirement in two dimensions that equal amounts of fiux must be reconnected 

\^ ' from each upstream region. The X-line drifts away from the conducting wall in all 

simulations with asymmetric outflow and into the strong magnetic field region during 
_^ ■ most of the simulations with asymmetric infiow. There is net plasma fiow across the 

^^ ■ X-line for both the inflow and outflow directions. The reconnection exhaust directed 

^_J . away from the obstructing wall is significantly faster than the exhaust directed towards 

it. The asymmetric inflow condition allows net vorticity in the rising outflow plasmoid 
which would appear as rolling motions about the flux rope axis. 
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INTRODUCTION 



Flux rope models of coronal mass ejections (CM Es) predict the formation of an elongated cur- 
rent sheet in the wake beh ind the rising plasmoid (e.g.. lSturrocklll966l : lHiravamalll974l : lKopp Sz P 



1976 



neuman 



Lin &: Forbesll200Cll ). Reconnection in these current sheets increases the flux contai ned within 



the ris ing plasmoid and reduces the amount of flux confining the plasmoid to low heights. iLin et al 



(J2004l ) predict that these current sheets contribute substantially to the mass budgets of CMEs such 
that the final parameters for CME evolution for interplanetary propagation are not set until the 
flux rope reaches sev eral solar radii. CME current sheets ma,y also play an important role in th e 



CME energy budget (jRakowski et al.ll2nn7l : iReeves et al.ll2nid : lMurphvll2nid : lMurphv et al.ll2nil 



In recent years, several features identified as current sheets have been observed during a nd after 



CMEs ( Ciaravella et al. 



2002; 



Ko et al. 



2003; 



Sui fc Holma^bood: 



Webb et al. 



2003; 



Lin et ah 2005,; 



Bemporad et al. 



2001 



a 



Ciaravella fc Raymond 



2008 



1 



Bemporad 



2008 



Lin et al. 



2OO9I: Vrsnak et al. 



200s 



2010 



Saint- 



Hilaire et al.ll2009 



Schettino et al 



201 



Liu et al 



2013; iReeves fc Golubl I2OIII : IPatsourakos fc VourlidasI I2OIII ) . These features usually appear 



20101 : ISavage et al.ll201C 



2011 



Landi et al. 



as long-lasting bright streaks when viewed edge-on in white light observations by the Large Angle 
and Spectrometric Coronagraph (LASCO, iBrueckner et al.l Il995l ) on the Solar and Heliospheric 
Observatory (SOHO). Several even ts observed by L ASCO were also observed by the Ultraviolet 



images occ asionally show large blobs (e.g, 



plasmoids (JFermo et al.ll2010l : lUzdenskv et al.ll2010l : 



Coronagraph Spectrometer fUVCS. iKohl et al.lll995) on S OHO in [Fe xviii] and Si xii. LASCO 



Lin et al. 



2OO5I) that are poss ibly the result of merging 



Loureiro et al.ll201lh . although the UVCS de- 



tection of strong C III emission in a few of these features has indicated that this low ionization state 
plasma was not processed through the current sheet. While these blobs have been interpreted as 
propagating at or near the Alfven velocity, Va, simulations show that large magnetic islands often 
propagate at a velocity ~2-4 times slower (e.g., IShen et al.ll201ll ). Frequently, these current sheet 



features are observed to drift o r tilt with time, including during the 'Cartwheel CME' (ISavage et al. 



2010 



Landi et al. 



201c 



2013) observed by the X-Ray Telescope (XRT) on Hinode on 2008 April 



9. The magnetic field topology and field strengths in and around CME current sheets are not well 
understood because of the lack of appropriate diagnostics. 

Magnetic reconnection with asymmetry in the inflow d irection occurs in the Ea rth's dayside 



magn etopause (jPhan Sz Paschmannlll996l:lKu Sz Sibec 



20041 ). la boratory plasina experime nts (JYamada et al. 



1997 



19971'). the Earth's magnetotail (^Oieroset et al 



Murphy fc Sovinedl2008), plasma tur - 



bulence ( Servidio et al.ll2009l . l20ld ). and during the merging of uneq ual flux ropes (LintonI 20061) . 



2008 



The scaling of asymmetric inflow reconnection has been investigated bv lCassak Sz Shavi (|2007l . 
20091 ) who find that the reconnection rate depends on a hybrid Alfven speed that is a function of 
the density and magnetic field strength in both upstream regions, given in dimensionless form by 



Vlk 



BlBr{Bl + Br) 



(1) 



PlBb^ + prBl 
where Bl and Br are the upstream magnetic field strengths and pi and pR are the upstream den 



sities (see also iBorovskv &: Hessd 120071 : iBirn et al.ll2008l . l20ld ). Simulations persistently show that 
the flow stagnation point and magnetic field null are separated by a s hort distaiice, and that the 
X-liiie usually drifts toward th e region with the stronger magnetic field (|Ugaill2000l : ICassak &: Shav 



20071 : iMurphv &: Sovinedl2008l ) . The current sheet drifting is less preval ent when a spatially local 



ized resistivity is imposed (e.g.. lBorovskv Sz Hessdl2007l : iBirn et al.ll2008l ). When there is a pressure 
gradient along the infiow direction during guide field reconnection, the X-line diamagnetically drifts 
along the outflo w direction. Reconnection is suppressed when the drift velocity is comparable to the 
Alfve n velocity (jRogers &: Zakharovlll995l : ISwisdak et al.ll2003l : |Phan et al.ll2010l : lBeidler &: Cassak 



2011 



CME current sheets form in a stratified atmosphere and thus are expected to be asymmetric 
along the outflow directio n as well. Asyrnmetri c outfiow reconnection occurs i n the solar atmosphere 



2006 



tails (Oka et al.l 



Lin et al. 



(IKoDP fc PneumanI 1 19761 : iLin fc Forbes 



200C 



20111 ') , and laboratory plasma e xperiments (lOno et al.lll993l. 11997 ; llnoinoto et al 



20081 : iMurphv fc Sovined l2008l : lOrav et al.ll20inl V iMurphv et al. 



Ciaravella &: Ravmondl 20081) . plan e tary magneto- 



(|20inl l derived 



scaling relations for a long and thin current sheet with asymmetric downstream pressure and 
found that the reconnection rate is not greatly affected unless outflow is blocked from both ends 
of the current sheet. Simulations of asymmetric outflow reconnection and X-line retreat show 
that the X-line and flow stagnation point are separated by a short distance and that most of 
the energy released during reconnection is directed towards the unobstructed current sheet exit 



(JRoussev et ahlboOll : ICalsgaard fc Roussevl[20o3 : loka et al.l[2008l : iReeves et allboiol : lMurphvl[201ol : 



Shen et al 



2011 



Asymmetry a nd motion of magnetic nulls have recently been considered in three- 



dime i isional configurations (e.g. JAI-Hachami &: Pontinll2010l : iGalsgaard fc: Pontinll201ll : iGray et al 



2010; Lukin fc Linton 2011 



In this work we simulate magnetic reconnection that is asymmetric in both the inflow and 
outflow directions and consider the effects that these asymmetries may have on solar flare and 
CME current sheets. Li Section [2l we discuss the numerical method used by the NIMROD code 
and the simulation setup. In Section [3l we present the simulation results including the magnetic 
structure, flow pattern, X-line dynamics, vorticity in the outflow plasmoid, and the morphology of 
the post-flare loops. In Section [4l we derive an analytic solution for the potential fleld structure 
of asymmetric post-flare loops and discuss the connection with the simulation results. In Section 
[5l we discuss expected observational signatures of line-tied asymmetric reconnection during solar 
eruptions. Section [6] contains a summary and a discussion of our results. 



2. NUMERICAL METHOD AND SIMULATION SETUP 



The NIMROD code (jSovinec et al.ll2004l . l2005l . l2010l ) solves the equations of extended magne- 
tohydrodynamics (MHD) using a finite element formulation for the poloidal plane, and for three- 
dimensional simulations, a finite Fourier series expansion for the out-of-plane direction. In di- 
mensionless form, the equations solved for the two-dimensional simulations reported in this paper 



are 



,'dv 



_ = -V X (r?J - V X B) , 
J = V xB, 

= JxB-Vp + V- pi/VV, 



P 



7- 



.vf— ) 



-|V-V-V-q+Q, 



(2) 

(3) 
(4) 

(5) 
(6) 



where B is the magnetic field, J is the current density, V is the bulk plasma velocity, p is the plasma 
pressure, p is the density, r] is the resistivity, z^ is the kinematic viscosity, D is an artificial number 
density diffusivity, T is the temperature, 7 = 5/3 is the ratio of specific heats, and Q includes Ohmic 



and viscous heating, 
and is given by q = 



The heat flux vector inc^ 



X||bb + x± (I-bb 



udes the effects of anisotropic thermal conduction 
• VT, where b is a unit vector in the direction of 



the magnetic field. The normalizations are given by: -Bo, Po, Lq, to, Vao = Bq/^pqp = Lq/Iq, 
po = Bq/po = poV^f), Jo = Bq/pqLo, and ryo/^o = j^n = Yn = Dn ^ L l/tn. Divergence cleaning 
is used to prevent the buildup of divergence error (jSovinec et al.ll2004l ). A small number density 
diffusivity {D <^ r]) is included to ensure that the number density profile remains sufficiently smooth 
in regions of sharp gradients. Except as noted, all velocities are in the simulation reference frame. 
A dditional details for the numerical method, normalizations, and simulation setup are presented 
bv lMurphvl(|2O10l ). 



Rather than modeling a solar eruption in detail, we choose an idealized initial condition that 
allows asymmetric reconnection to commence. The initial conditions are of a perturbed, modified 
Harris sh eet that allows a s ymme tric upstream magnetic field strengths, densities, and plasma pres- 



sure see 



Birn et al.ll2008l . I2OIOI ). We define x as the infiow direction, y as the outflow direction, 



and z as the out-of-plane direction. The initial equilibrium is given by 



B(x) = Bro 



tanh 



b]+b 



l + b 



y + Bzoz 



p{x) = \{1-B'^)+I3ro^ 
p{x) = po[l + {fo-l){l-X')], 



A = 



1 



1 + tanh 



X 

To 



(7) 

(8) 
(9) 

(10) 



where 6q is the initial current sheet thickness, b is the magnetic asymmetry parameter, Bzo is the 
initial guide field (equal to zero except in case F where B^o = 4), and Pro = pro/ (-B^q/2). Note 
that /3ro does not include the guide field contribution. The subscripts 'L' and 'i?' refer to the 
asymptotic magnitudes of quantities for x < and a; > 0, respectively, and the subscript '0' refers 





Bp = 


V X (Apz) , 


Ap = 


—Bphexp 


\h) \ h ) 



to the values of quantities at the beginning of each simulation. The initial ratios for the upstream 
densities and magnetic fields are given by 

/o = ^, (11) 

PRO 

^ Bu^^l-b 

5/jo 1 + fe ^ ^ 

where < 6 < 1, Blo > 0, and Brq > 0. These ratios at t = will in general differ somewhat from 
the ratios during the course of each simulation. 

The initial magnetic perturbation is of the form 

(13) 

where 

(14) 

Here, Bp is the strength of the perturbation and h governs the perturbation width. The perturbation 
is centered about {x,y) = (0, A). 

The simulations we report on are presented in Tabled) The controls for our study are case A, 
with symmetric inflow and asymmetric outflow, and case H, with asymmetric inflow and symmetric 
outflow. The eight simulations in Table [1] are chosen to test different upstream magnetic field 
asymmetries (compare cases A-D), different resistivities (cases C and E), the inclusion of a guide 
field (cases C and F), and different upstream densities (cases A and G). Case C is a control for 
comparisons with cases E, F, and H which all have Rq = 0.25. 

The simulation parameters are as follows. For the initial equilibrium we use Bq = 1, 6q = 0.1, 
and /3ro = 0.18. The initial perturbation is given by Bp = 0.1, A = 1, and h = 0.5, except for 
case H where A = 0. The diffusivities are given hy rj = u = 10^'^, D = 10~^, x± = 10~^, and 
X|| = 10~^. The Lundquist number must be considered carefully because the Alfven speeds differ 
in each upstream region. For our comparisons we use a hybrid Lundquist number based on the 
hybrid Alfven speed presented in Equation ([T|), 

^h^ , (15) 

V 

where L is the characteristic length scale of the current sheet. Mesh packing is used to concentrate 
resolution in regions of strong gradients. However, because the current sheets drift in most simula- 
tions, high resolution is needed over a larger part of the domain than in simulations with symmetric 
inflow and/or outflow. We use seventh order flnite element basis functions in all simulations. There 
are m^ and rriy flnite elements along the inflow and outflow directions, respectively. The size of 
the computational domain for cases A-G is given by — 7 < x < 7 and < y < 30. We assume 
conducting wall outer boundary conditions in all directions for cases A-G; consequently, late in 



time there is some influence from line-tying along the upper boundary at y = 30. However, we 
concentrate on dynamics far from the unphysical upper boundary (y < 18). Case H differs in that 
— 15 < y < 15 and there are periodic boundary conditions in the y-direction. While we include 
anisotropic thermal conduction, we neglect radiative losses, coronal heating, and vertical stratifi- 
cation of the atmosphere. The cooling time scale for CME current sheet plasma (T ~ 5 x 10^ K 
and n ~ 5 X 10^ cm^^) is about half a day. This is an order of magi iitude longer than a dynamical 



time scale so we are justified in neglecting radiative losses (see also llmada et al. 



2011 



It is appropriate to remark upon our choice of a uniform, explicitly defined resistivity. Occa- 
sionally, prior simulation s of asymmetric inflow r econnection have used a spatially localized resis- 



tivity enhancement (e.g.. lBorovskv fc Hessd 120071 ). The intended effect is to constrain the position 



of the X-line. However, this constraint is artiflcial and likely changes the internal structure of 
the reconnection region and the sca l ing b ehavior. We speculate that this is the reason why the 
scaling found by iBorovskv fc Hessd (J2007l ) differs slightly from the analytic prediction made by 
Cassak &: Shayi ()2007l ). Some simulations use resistivity as a function of current density or other 
plasma paramet ers so that the position of the X-line is not artificially constrained (see, for example, 
Birn et al.lll996l ). Alternatively, some simulations of reconnection have used numerical resistivity 
inherent in the discretization as the sole means of breaking the frozen-in condition, rather than an 
explicitly defined resistivity or field line breaking mechanism incl uded in the physical model (e.g., 



Laitinen et al. 



2005 



Ouellette et al. 



2010 



Edmondson et al.ll2010l ). Such models are effectively us- 



ing an ideal MHD algorithm. However, when there is significant numerical reconnection in an ideal 
MHD simulation, the results are by defini tion not converged. T heref ore, while scali i ig mo dels for 
asymmetric reconnection such as those by ICassak fc: Shayl (120071 ) and iMurphy et al.l (120101 ) do not 
explicitly depend on the dissipation mechanism, the lack of convergence makes ideal MHD simu- 
lations poor tests of asymmetric reconnection models. Our choice of a uniform, explicitly defined 
resistivity avoids artificially constraining the position of the X-line and allows convergence in the 
numerical simulations. 



3. SIMULATION RESULTS 



In this section, we describe the principal results of our simulations of line-tied asymmetric 
reconnection. When making comparisons between the simulations presented in Table [T] it is im- 
portant to note that the scaling for cases A-D is performed by keeping Brq constant and reducing 
Blq. Consequently, the total magnetic energy available to be reconnected and the hybrid Alfven 
speed given in Equation ([1]) decrease as the simulations become more asymmetric. Case G is best 
compared directly to case A, and cases E, F, and H should be compared to case C. 



3.1. Principal Features 

The general features of our simulations of line-tied asymmetric reconnection are presented in 
Figure dl which shows case C with an initial upstream magnetic field ratio of Rq = 0.25. Magnetic 
flux contours in Figure HJa) show that field lines in the strong magnetic field region are much 
less bent than field lines in the weak magnetic field region. The reconnected loops near the lower 
boundary and the outflow plasmoid both preferentially develop into the weak magnetic field region. 
The separatrices are traced by regions of strong out-of-plane current density in Figure mb), but the 
portions that bound the strong field regions have much stronger current density than the portions 
that form the boundary between the outflow and weak field regions. The plasma pressure buildup 
is almost entirely contained on the weak field side of the current sheet [Figure [T]^e)]. 

The upward outflow velocity, shown in Figure [TJd), is substantially faster than the downward 
outflow velocity because of the obstructing wall along y = 0. This behavior is consistent with pre- 



vious simulations of asymmetric outflow reconnection (IRoussev et al. 



2001 



Galsgaard &: Roussev 



2002; lOka et al.ll2008l : iMurphv fc Sovinedl2008l : lMurphvll2010l : iReeves et al.ll2010l l. This difference 
in velocities occurs for two principal reasons. First, outflow towards y = is obstructed by the 
buildup of plasma and magnetic pressure. Second, the X-line is located near the lower exit of the 
current sheet so that the tension force below the X-line directed towards y = is much weaker 
than the tension force above the X-line directed away from y = 0. 

In Figure Hl^d), there is a stream of plasma flow to the right of the X-line but extending 
downward into the post-flare loop structure and tracing the separatrix. The magnitude of the 
plasma flow is small (Vy ~ 0.02 compared to Vaho = 0.25), but is positive in a region where negative 
Vy is expected. Magnetic tension in this region is pulling plasma downward, but is countered by a 
comparable contribution by the vertical plasma pressure gradient and a modest contribution by the 
vertical magnetic pressure gradient pushing plasma upward. The current density is strong so that 
resistive diffusion acts to make the field lines more potential even against this plasma flow. Such 
flows might be observable, but three-dimensional geometry may make identification ambiguous, 
and it is not clear that they will occur when /3 ^ 1. 



3.2. Internal Structure 



The X-line (the line about which the magnetic field has a hyperbolic or X-line topology) is 
located at x„ = (x„, y„). This position is given as a function of time for several simulations in Figure 
[2j Figure El^a) shows the X-line position along the infiow direction. In most previous simulations of 
asym metric infio w reconnection, the X-line drifts into the region with the stronger magnetic field 
fe.g..lUgaill2000l ). This drifting is tied locally to spatial derivatives in the out-of-plane electric field 
(JMurphvll2010l ). In particular, the X-line drifts along the inflow direction towards decreasing E^. 
For many but not all of the cases we observe the X-line drifting towards the strong field upstream 
region for most of the simulation. In case D with Rq = 0.125, the X-line drifts into the weak field 



region before drifting into the strong field region. Case E is identical to case C except for having 
three times the resistivity, and shows a significant drift of the X-line into the weak magnetic field 
region. Comparing cases A and G suggests that the inclusion of a density asymmetry does not lead 
to a significant drift along the inflow direction. The rate of X-line drifting along the inflow direction 
is a function of the magnetic field asymmetry, the resistivity, the inclusion of a guide field, and the 
simulation setup. However, the dependences on each of these parameters are not straightforward. 

The X-line position retreats from the obstruction at y = for all asymmetric outflow cases 
[Figure [D^b)]. The rate of X-line retreat is comparable for cases A and B (with Rq = 1 and Rq = 0.5, 
respectively) but decreases as the initial conditions become more asymmetric. By comparing cases 
C and E we see that increasing the resistivity actually slows down X-line retreat, in part because 
it results in smoother gradients in E^- From cases A and G we see that increasing /o from 1 to 4 
decreases the rate of X-line retreat slightly. Comparing cases C and F shows that including a guide 
field can lead to much faster X-line retreat. In contrast to all other asymmetric outflow simulations, 
the X-line is located near the top exit of the current sheet in case F. 

Previous simulations of reconnection with either asymmetric inflow or asymmetric outflow have 
generally shown a separation between the flow stagnation point and the principal X-line when a uni 



forrn resistivity is used fe.g.. ICassak Sz Shavll2007l . 120081 . 



2008 



Murphy 



2010 



Shen et al. 



2011 



200S 



Oka et al.ll2008l : lMurphv fc Sovinec 



Contrary to our expectations, there is in general no flow 
stagnation point associated with the X-line in the simulation reference frame for cases with both 
asymmetric upstream magnetic fields and asymmetric outflow. There are contours where V^ = 
and Vy = 0, but these contours do not intersect with each other in the current sheet for most of 
the simulations with asymmetry in both the inflow and outflow directions. There still is a reversal 
of the outflow component of velocity near but not colocated with the X-line. 

Figure El shows both the plasma flow velocity at the X-line, V„ = {Vx {xn,yn) ,Vy {xn,yn)), 
and the rate of change in position of the X-line, x„ = {xn,yn), for case C with Rq = 0.25. The 
plasma flow velocity at the X-line differs greatly from the X-line drift velocity along both the inflow 
and outflow directions, indicating a signiflcant departure from the frozen-in condition. For the 
inflow direction, Vx{xn,yn) is of the same sign as Xn, but Vx{xn,yn) is signiflcantly greater. For 
the outflow direction, y„ remains positive but Vy{xn,yn) becomes negative so that the X-line is 
retreating against the flow of the plasma. The X-line is able to diffuse a gainst strong pl asma flow 
by diffusion of B x along the inflo w direction, as described in detail by iMurphvl (J2010l ) (see also 



Siscoe et al. 



2002 



Oka et al.ll2008l l. 



A slice along the inflow direction for case C at the height of the X-line {y = yn) is given in 
Figure S] for t = 100. The maximum in J^ is located on the strong fleld side of the X-line. This 
result that the X-line is on the weak fleld side of the current sheet is consistent with the simulations 



of lCassak &: Shavl (|2007l ) . The X-line is located near a local maximum in plasma pressure along the 
inflow direction. In the simulation reference frame, Vx ^ away from the diffusion region on the 
strong fleld side of the current sheet. This can be interpreted as most of the inflow coming in from 



the weak field side, with the X-hne and current sheet drifting into the strong field upstream region. 
The outflow component of velocity, Vy, has an interesting bipolar signature. While the X-line is 
moving in the positive y direction, the X-line is located near where Vy is most negative. The inflow 
component of the magnetic field, Bx, is negative throughout the slice except at the X-line where 
it is zero. This profile allows negative B^ to dif fuse inward so that at later times the X-line is 



located at higher heights (compare to Figure 5 of lMurphyll201Cl ). The reconnecting component of 
the magnetic field. By, shows that the initial ratio of Rq = 0.25 is somewhat less extreme than the 
ratio oi R ^ 0.17 at t = 100. 



3.3. Reconnection Rate 

The reconnection rate, defined as the out-of-plane component of the electric field at the X-line, 
is shown in Figure [5] for each of the cases described in Table [TJ Comparing cases A-D shows that 
decreasing the magnetic asymmetry factor Rq leads to a corresp onding decrease i i i the reconnection 



rate. This is qualitatively consistent with the scaling derived bv lCassak &: Shavi (|2007l ). Increasing 
the density asymmetry factor /q decreases the reconnection rate modestly from case A to case G. 
The reconnection rate in case C is somewhat quicker than case E, even though case E is three 
times more resistive. The inclusion of a guide field slightly increases the reconnection rate in case 
F compared to case C. The reconnection rate is modestly quicker with asymmetric infiow and 
outfiow (case C) than in an otherw i se equ ivalent simulation with asymmetric inflow but symmetric 



outflow. As described bv iMurphvl (J2010l ). this occurs because the current sheet can only increase 



in length along one outflow direction. 



3.4. The Outflow Plasmoid 

The asymmetric inflow condition allows net vorticity in the outflow plasmoid. In Figure [6] we 
show velocity vectors in the reference frame of the 0-point for case C at t = 100. The counter- 
clockwise flow pattern is largely due to the reconnection outflow jet impacting the rising plasmoid 
at an angle; consequently, the outflow jet increases the net vorticity in the plasmoid. This circu- 
lation should be a generic feature of asymmetric inflow reconnection but not asymmetric outflow 
reconnection. This interpretation should be considered qualitatively rather than quantitatively 
because the presence of closed field lines is partially due to reconnection near the outer, artificial 
line-tied boundary at y = 30. However, these results suggest that the orientation of reconnection 
outfiow jets relative to the fiux rope can lead to an analogous circulation pattern when the outfiow 
jet impacts the fiux rope obliquely rather than directly at its base. 

In contrast to the X-point, the frozen-in condition near the 0-point is approximately met. 
For example, at t = 100 in case C, the plasma fiow at the 0-point is approximately equal to 
the velocity of the 0-point: {Vx{xo,yo),Vy{xo,yo)) = (—0.0076,0.1211) compared to {xcijo) = 



10 



(—0.0065,0.1238), where Xq = {xo,yo) is the position of the 0-point. This indicates that the O- 
point is being primarily advected by the bulk plasma flow and that diffusive flow across the 0-point 
is not signif i cant. Plasma flow across an 0-point can occur in a process similar to that described 



by iMurphvl (|20ld ) for X-line retreat. When the 0-point in a magnetic island is displaced towards 
one particular direction, resistive diffusion will in general act to change the position of the 0-point 
to be closer to the center of the island when the island is not significantly distorted. 



3.5. Morphology of the Post-Flare Loops 

The most easily observable difference between the symmetric and asymmetric cases is the 
structure of the post-flare loops. In particular, the asymmetric loops are skewed when compared 
to the symmetric case, and take an asymmetric candle flame shape. While in the symmetric case 
the tops of each loop are all along x = 0, the loop-tops in the asymmetric case have their apexes at 
different locations along the inflow direction (Figure [Tj) . At low heights where the field has relaxed 
to a near-potential state, the loop-tops are displaced toward the low magnetic field side. At greater 
heights, the loop-tops become located closer to the current sheet demarcating the low and high 
magnetic field regions. 



4. AN ANALYTIC SOLUTION FOR THE POTENTIAL FIELD STRUCTURE 

OF ASYMMETRIC POST-FLARE LOOPS 

In Section [331 and Figure [7] we show that the location of the post-flare loop apexes in each case 
are a function of height. In the current section, we present an analytic solution for the potential 
field structure of post-fiare loops in the present configuration. This solution provides insight into 
the observational signatures very near the magnetic field reversal along the lower boundary, but 
should be treated as a limiting case for post-fiare loops late in time that have been able to relax to 
a near-potential state. 

We consider the domain < x < vr and y > 0. The boundary condition along y = is given 

by 

By{x,0) = < , (16) 

nji : a < X < TT 



where the constants Bl and Bji are positive and the location of the field reversal is given by 

_ ttBr 
""' Bl + Br 



(17) 



Consequently, there is no net magnetic fiux from the lower boundary over < a; < vr. Because the 
outer boundary is artificial, we must consider regions close to the magnetic field reversal to avoid 
these effects. However, this analysis also assumes that the magnetic field reversal length scale is 
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smaller than the region being investigated. Because we are finding a potential field solution we also 
require that /3 ^ 1, or more generally, that plasma pressure gradient forces are small. 



The potential magnetic field is given by 

B = -Vtp, 
where the scalar potential ^ is governed by Laplace's equation, 

0. 
The unique solution to Laplace's equation appropriate to our boundary conditions is 



vV 



(18) 



(19) 



^ = E 



2 [Bl + B 



ra=l 



vr 



n^ 



* sin na cos nx e "^. 



Equation (j20p may be written in the form of a vector potential B 
expression for the magnetic flux, 



A., 



2 [Bl + Br 



n=l 



vr 



n^ 



sm na sm nx e 



-ny 



(20) 
V X {Azz) to provide an 

(21) 



On large scales, the solution depends on the outer boundary. On scales much smaller than the outer 
boundary, the solution becomes scale- free as we approach (x, y) — )• (a, 0). Because in general Fourier 
series expansions will be truncated at some A^max, we must also consider scales Ax S> vr/A^max- Using 



Euler's formula and the identity In (1 
are given in closed form by 



J2'^=i Q^'/n, the components of the magnetic fleld 



B-r 



Bn 



Bl + Br 



In 



27r 

Bl + Br 

2tt 



(1 



J,a+ix—y 



)(1 



-ia—ix—y\ 



(1 



In 



j,a—ix—y 



(1 



-ia+ix—y 



1 



^ia+ix—y 



) (i - e^-i^-y) 



[l — Q-ia+ix-y\ n 



-V) 



(22) 
(23) 



Solutions for four different magnetic field ratios are presented in Figure [8] very close to the magnetic 
field reversal along the lower boundary. These solutions approximate the magnetic field structure 
of post-flare loops after these loops have had time to relax. These loops have a similar appearance 
to the reconnected loops very near the fleld reversal that are shown in Figure [71 

Next we derive the loop-top positions along the inflow direction as a function of height. By 
setting By = and using Euler's formula and trigonometric identities, we can obtain an equation 
for the height of loop-tops as a function of x. The loop-top positions are given by 



Vlt{x) = In (cos a sec x) . 



(24) 



For our purposes this expression is most useful near the fleld reversal along the lower bound- 
ary at (x, y) = (a, 0) where the solution is insensitive to the outer boundary. We deflne 6 = 
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arccot (dy^j'/dx) as the (clockwise) angle with respect to vertical for the loop-top positions. Be- 
cause dyir/dx = tanx and R = Bl/Br, evaluating Equation (p^ at x = a yields 



vr /R- 1\ , , 

(25) 



2 V^+1 

Equation (j25p provides an upper limit on the distortion of the asymmetric post-flare loops as a 
function of asymmetry. The angle from vertical for the loop-top positions is shown in both Figures [7] 
and El In Figure [7] we see that Equation (I25p reliably approximates the angle from vertical that the 
loop-top positions take very near the field reversal. Equation (125p overestimates the angle slightly 
in the simulations because of the finite width of the magnetic reversal along the lower boundary 
((^0 = 0.1) and because the field lines are not fully potential due in part to a finite pressure gradient. 



5. OBSERVATIONAL SIGNATURES OF ASYMMETRIC RECONNECTION 

In this section, we describe the observational signatures predicted by our simulations of line- 
tied asymmetric reconnection during solar eruptions. 



5.1. The Location of the X-line and Flow^ Reversal 

The location of the principal X-line is important in flare/CME current sheets because it helps 
determine the partition of outflow energy towards and away from the Sun. These simulations 
suggest that, at least during reconnection with a weak guide field, the principal X-line will be 
located near the lower base of fl are/CME current sheets, not too far above the post-flare loops (e.g., 



Murphvll20ld : IShen et al.ll201ll ). The principal X-line and flow reversal are probably separated by 
a distance that is shorter than the observational errors and systematic uncertainties. Finding the 
location of the flow reversal in a current sheet viewed nearly edge-on is difficult but can be done 



by tra cking the motions of current sheet blobs. This task has been accomplished bv lSavage et al 



(|201Cl ). who find that the flow reversal is at a height of just ~0.25 solar radii above the limb in the 
Cartwheel CME current sheet. For comparison, the post-flare loops have a height between 0.1 and 
0.2 solar radii above the limb and the current sheet extends several solar radii outward into the 
LAS CO held of view. Our simulations also predict that the X-line is located on the weak magnetic 
field side of the current sheet (see Figure [I]); however, we anticipate that this signature is beyond 
our current observational capabilities. 



5.2. Post-Flare Loops and Loop Footpoints 

A signature of line-tied asymmetric reconnection is the distortion or skewing of post-flare loops 
as shown in Figure [7] and described in Section [3.51 The loop apexes are not immediately above each 
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other. Rather, the apex positions along the inflow direction are a function of height. The post-flare 
loops develop a characteristic skewed candle flame shape. This signature should be apparent in 
Ha, EUV, and X-ray observations of post-flare loops during line-tied asymmetric reconnection. A 
candidate event with a cle ar skewed cand l e flam e sh a.pe is the 1992 February 21 flare on the east 
limb that was analyzed bv iTsuneta et al.l (jl992l ) and iTsunetal (|l996l ). However, projection effects 
associated with a complicated three-dimensional geometry m ight also lead to a skew ed candle flame 
structure during some events (see, for example. Figure 15 of lForbes &: Actonlll996l ). 



Solar flares characteristically show hard X-ray (HXR) emission at the footpoints of newly re- 
connected field lines in response to energetic particles and the thermal conduction front impacting 
the chromosphere. The standard model of solar fiares predicts that the footpoints move away from 
the neutral line so th at the instantaneous location of the footpoints i s given by the amount of 



reconnected fiux (e.g., iLin et al. 



1995 



Forbes &: Acton 



1996: 



Lin 



observed during many flares (e.g., lAsai et al.l l2004l : iKrucker et al 
although more coniplicated motions are possible (e . g^ 



Grigis &: Benz 



2005 



Somov et al. 



2005 



Ji et al. 



2006 



2004l\ . This behavior h a s bee n 



2005; 



Bogachev et al. 



Su et al. 



W. 



Yang et al. 



200,51 : 



2009, 2011 



Sakao et al. 



199S 



20071 ) . Because of the requirement 



that equal amounts of flux be reconnected from each upstream region in two-dimensional simula- 
tions, we predict that the velocity of the footpoint in the strong fleld region will be slower than the 
velocity of the footpoint in the weak fleld region. For example, if the magnetic field in one upstream 
region has twice the strength of the other region, then the footpoints on the strong magnetic field 
side will move half as quickly as the footpoints on the weak magnetic field side. By combining 
observations during the 2003 October 29 XIO flare by the Reuvan Rama ty High Energy Spec tro- 
scopic Imager (RHESSI) an d the Michelson Doppler Interferometer (MDI, IScherrer et al.lll995l ) on 
SOHO, iKrucker et al.l (|2005l ) showed that the magnetic fleld n ear the slower moving footpoint was 
generally stronger than in the faster moving footpoint (see also ISvestkalll976l ) . While this is consis- 
tent with the predictions from our models, theory needs to take into account the three-dimensional 
nature of magnetic reconnection and the patchy distribution of magnetic flux in the photosphere. 

Asymmetry in the post-flare loop structure will affect the relative intensities from the emis- 
sions at each footpoint. HXR emission is largely determined by the transport of energetic particles 
from above the loop-top into the chromosphere. Energetic particles on the strong magnetic field 
side will be more likely to be reflected because magnetic mirroring is more effective. Therefore, 
energetic particles will be more likely to enter the lower solar atmosphere on the weak magnetic 
field side. Consequently, the footpoint on the weak magnetic field si de is expected to yield stronger 
HXR emission and chromospheric evaporat io n due to bearn . heatin g (JMelrose fc Whitelll979l : I Sakao 
E994I : iMehose k Whitel Il98ll : iHuangI boOTJl . iKundu et al.l (|l995l ) present two flares observed by 
the Nobeyama rad ioheliograph ([Nakaiima et al.lll994l ) and the Hard X-ray Telescope on Yohkoh 
( Kosugi et al.lll99ll ) in support of this scenario. The footpoints with weak HXR emission (indicat- 
ing fewer particles escaping from the trap) had stronger gyrosynchrotron e mission (suggest ive of a 



stronger magnetic fleld). While this behavior does not happen in all flares (iGoff et al. 



tional sources of asymmetry include an asymmetric initial pitch angle distribution (e.g.. 



2004), addi- 



Huang et al. 
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20ld ) , differences in the column density in each fo qtpoint (see, however, ISaint-Hilaire et al.1 12008| ) 
and directionahty in the accelerating electric field (|Hamilton et al.ll2005l : iLi &: Linll2012l ). 



5.3. Reconnection Inflow Velocities 



Several works have report ed observations of reconnection inflow velocities associated with 
flare/CME current sheets (e.g. , lYokovama et al.l I2OOII : iLin et al.l l2005l : iNarukage fc Shibatal 120061 : 
Takasao et al.ll2012l ). During asymmetric inflow reconnection, the inflow velocities on either side of 
the current sheet are predicted to differ. In a steady-state, one would expect from flux conservation 
that the out-of-plane electric field will be constant, leading to the relation 



VlBl = VrBr. 



(26) 



This implies that the ratio of inflow velocities is inversely proportional to the ratio of upstream 
magnetic field strengths and that in principle measuring the ratio of infiow velocities on either side 
of the current sheet would directly provide the upstream magnetic field ratio. 

Figure U] shows that there are strong variations in the electric field in the simulation reference 
frame and that the inflow component of velocity is approximately zero in the strong upstream region. 
Thus the reconnection process in the simulations is not time-independent. Therefore, Equation (j26p 
should not be expected to give reliable estimates of the magnetic held asymmetry except by taking 
the velocities in the reference frame of the X-line and showing that the reconnection process is 
steady. However, despite the inapplicability of Equation (f26]l during time-dependent asymmetric 
reconnection, a systematic demonstration that the reconnection inflows differ on either side of the 
current sheet will provide suggestive evidence in future observations that the reconnection process 
is asymmetric. 



5.4. Drifting of the Current Sheet 

A commonly observed feature of CME current sheets is that they appear to drift or tilt with 
time. ffwode/XR T observations of t h e Cartwheel CM E show that the current sheet drifted at a 
rate of ~4° hr~^ (jSavage et al.ll2010l ). IKo et al.l (|2003l ) discussed observations on 2002 January 8 
of a current sheet with a drift of ~0.8° hr~^. These drift rates are signiflcantly greater than can be 
accounted for by solar rotation. 

There are seve ral possible explanations for the observed drifting of CME current sheets: (1) 
Savage et al.l (J2010l ) suggest that the Cartwheel CME current sheet is observed at an angle to the 
plane-of-sky and that the appearance of drifting is caused by different parts of the current sheet 
actively reconnecting at different times. (2) The tilting could be due to the drifting of the X-line 
and current sheet into the strong magnetic fleld region during line-tied asymmetric reconnection as 
discussed in this paper. The predicted velocity of less than a percent of the Alfven speed is within 
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observational constraints. However, the simulations show drifting but not the observed tilting of 
the current sheet. (3) The rising flux rope could pull the plasma sheet region along with it so 
that the drifting is caused by macroscopic behavior. This mechanism requires that the current 
sheet become more aligned with the direction of flux rope propagation. For the Cartwheel CME, 
however, this appears to not be the case. (4) There is a large-scale force imbalance between the 
two upstream regions so that the entire region surrounding the current sheet is pushed rapidly 
tow a.rds one directi on. Such a force imbalance could easily occur in the early stages of an ejection 
(see iKo et al.l 120031 ) . (5) The tilting is caused by relaxation in the post-eruption active region as 
the magnetic field configuration becomes more potential. 

Further numerical and observational tests are required to constrain which mechanisms lead to 
current sheet drifting during CMEs. The location of the CME current sheet relative to the direction 
of propagation of the flux rope may be of particular importance in determining how upflow from 
the current sheet impacts and influences the evolution of the rising plasmoid. However, it is not yet 
known if CME current sheets are an important component of the energy and momentum budgets 
of CMEs. 



5.5. Circulation Within the Rising Flux Rope 

In Section [3.41 we show that the asymmetric inflow condition allows net vorticity to develop in 
the outflow plasmoid. This circulation pattern develops because the outflow jet impacts the rising 
flux rope obliquely rather than directly at its base. When the current sheet and direction of flux rope 
propagation are misaligned, this leads to the possibility that the reconnection outflow jet torques 
the rising plasmoid. Vortex motions like those seen in Figure [6] are not unique to configurations with 
a line-tied lower boundary, but rather should be a generic feature of reconnection with asymmetric 
upstream magnetic fields such as at the dayside magnetopause. 



Martini (J2003l ) , iPanasenco &: Martini (J2008l ) , and iPanasenco et al.l (|201ll ) report observations 



of several CMEs that display a rolling motion about the axis of the erupting prominence. We 
hypothesize that this rolling motion is induced by an offset between the CME current sheet and 
the rising flux rope during some events. However, it is unknown if the kinetic energy released by 
reconnection is enough to drive this circulation. 

There are other candidate mechanisms for the development of apparent circulation and vorticity 
in the rising flux rope when viewed in cross-section. A similar flow pattern could develop by un- 
twisting of the magnetic field in the rising flux rope during expansion and relaxation. Alternatively, 
the external magnetic field coul d deflect the flux r ope's outward motion, thus leading to apparent 
rolling behavior. For example, IPanasenco et al.l (120111 ) note the presence of coronal holes near 
CMEs that display rolling motion. This rolling behavi or about the flux rope axis is in contrast to 
rotat i on about the direction of CME propagation (e.g.. lLynch et al.ll2009l : iPatsourakos &: Vourlidas 
201 ll : [Thompson et al.ll2012l ^. 
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5.6. Ultraviolet Spectroscopy of the Inflow Regions and Current Sheet 

If there is a strong asymmetry in the plasma conditions on the two sides of the current sheet, 
an asymmet ry in the e missio n hne intensities might be expected as well. An example is seen in 



Figure 12 of iKo et al.1 (|2003l ). where the [Fe XVIII] and [Ca XIV] lines are sharply peaked at the 
same position along the slit, while the Si XII emission is shifted about 70" (one bin along the 
entrance slit, or 50 Mm) toward the north, and [Fe XII] is shifted still farther north by another 70". 
Lower temperature lines such as [Fe X] are entirely absent from the immediate area of the current 
sheet, but a definite [Fe X] feature is seen about 140" to the south. 

There is a clear a symmetry in the temperature between the north and south sides of the current 



sheet in Figure 12 of iKo et al.l (|2003l ). It is hard to estimate the temperatures accurately because 
of the ambiguity caused by the uncertainty in the foreground and background contributions to 
the line intensities. Roughly speaking, the [Fe X] and [Fe XII] intensities on the southern side are 
comparable, indicating a temperature logT ~ 5.95, while on the northern side Si XII is several 
times stronger than [Fe XII], suggesting logT ~ 6.15. 

The density contrast is more difficult to estimate. Based on those temperatures, the emission 
measure on the northern side is four times that to the south. Assuming that the depth along the 
line of sight is the same on both sides of the current sheet, the density is twice as high to the north. 
Another density diagnostic is the Ly/3 to Lya ratio. Since Ly/3 contains comparable contributions 
from collisionally excited (oc n^) and radiative scattering (oc n) components, while Lya is almost 
entirely formed by radiative scattering, the ratio is a reliable diagnostic. Unfortunately, the sepa- 
ration of the emission from the region of the current sheet from broad background emission is even 
more difficult than it is for the Fe and Si lines discussed above. We estimate that the Ly/3 to Lya 
ratios are 0.01 and 0.006 on the southern and northern sides, respectively, which would imply that 
the density is ~1.5 times higher in the south than in the north. Each of these density diagnostics 
is accurate to within about a factor of two, suggesting that the northern and southern densities are 
comparable to within uncertainties. Overall, we conclude that the plasma pressure is somewhat 
higher in the north than the south. It would be dangerous to conclude that the magnetic pressure 
is higher in the south than the north, however, since the region is probably not in equilibrium. 

Asymmetry in the upstream temperatures might also affect the charge state distribution in the 
current sheet plasma and the rising flux rope. For example, the two upstream regions could start in 
ionization equilibrium at T^ = 1 and Tr, = 2 MK. Plasma entering the current sheet will be heated 
very quickly to ~5-10 MK. The ionization time scales typically range from 10-1000 s and are a 
function of both temperature and density. This is comparable to a dynamical time scale and thus w e 



2011 



expect the plasma leaving the current sheet to be underionized (jKo et al.ll2010l : llmada et al.l 
Because the current sheet plasma contains contributions from upstream plasma with different initial 
temperatures, this may result in a broader charge state distribution than would be expected if there 
was just one starting temperature. However, the difference may be small because the time scales to 
ionize low charge state species will be short in very hot plasma. Differences in the inflow velocities. 
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densities, and temperatures might also lead to observable differences in the non-thermal line widths 
between the two upstream regions. 



6. DISCUSSION AND CONCLUSIONS 



In this article, we present resistive MHD simulations of line-tied asymmetric reconnection in 
the context of solar flare and CME current sheets. There is asymmetry along the inflow direction 
because of different upstream densities and magnetic fleld strengths, and asymmetry along the 
outflow direction because the initial perturbation is placed near a conducting wall. The simulations 
are used to understand the basic physics of asymmetric reconnection and predict observational 
signatures of asymmetric reconnection in the solar atmosphere. This approach provides a unifled 
picture of how asymmetry in the reconnection process impacts flare emissions and CME evolution. 

As in previous simulations of asymmetric outflow reconnection, the X-line is generally lo- 
cated near the lower base of the current sheet. Consequently, the outflow velocity in the un- 
obstructed direction is significantly faster than the outflow velocity in the obstructed direction. 
The slow downflows occur in simulations for two reasons. First, the downflows impact a re- 
gion of strong magnetic and plasma pressure, thus slowing them down. Second, the principal 
X-line is often located near the lower base of the current sheet so that the downward tension 
force is much weaker compared to the upward tension force. Downflowing loops in flare/CME 
current sheets are frequently observed to propagate at velocities several times slower than the 



Savage fc: McKenzie 



2011 



Warren et al. 



Alfven speed (IMcKenzie fc Hudson! Il999l: lAsai et al.1 12004: ISheelev et al.ll2004l : iReeves et al. 



2011 



McKenzi 



2011 



2008 



), in contrast to symmetric models of 
reconnection that predict bidirectional Alfvenic jets. Simulations such as those presented in this 
pap er suggest that the outflow asymmetry is responsible for this difference in outflow speeds (see 



also Roussev et al.ll200ll : ICalsgaard fc Roussevll2002l : lReeves et al.ll20ld:lMurphvll2010l : IShen et al 



20 111 ). Observational analysis techniques such as those by lSavage et al.l (l2010l ) that constrain the 
position of the reconnection flow reversal are important for resolving the problem of slow down- 
flows. ^We_do_jiote_that there are observations of fast downflows during some reconnection events 



(e.g., Ilnnes et al.ll2003l ). 



In simulations with asymmetric upstream magnetic flelds, we flnd that the post-flare loops 
are skewed when compared to cases with symmetric inflow. The loop-tops are not directly above 
each other; rather, the positions of the loop-tops are a function of height. The observed shape 
of these loops is reminiscent of a candle flame. The structure of near-potential loops close to 
the field reversal along the lower boundary are approximately given by our analytic solution in 
Section [H Cuspy p ost-flare loops with an apparent candle flame structure have been observe d 
during many events (JTsuneta et al.lll992l : lTsunetalll996l : iForbes &: Actonlll996l : iReeves et al.ll2008l ). 
Observations showing such behavior would allow us to diagnose or place limits on the asymmetry 
in the upstream magnetic fields. However, 3D geometry and projection effects would need to be 
considered carefully because they could also contribute to a skewed post-flare loop appearance (see 
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Forbes fc ActoiJll996l ). 



The morphological features associated with line-tied asymmetric reconnection grow preferen- 
tially into the weak magnetic field upstream region. Unlike in cases with symmetric inflow, there 
is net vorticity in the outflow plasmoid because the reconnection outflow jet impacts it obliquely. 
This result suggests that CME current sheets could drive a circulation pattern in erupting flux 
ropes when viewed in cross-section for events where the current sheet is offset from the direction 
of flux rope motion. 

The X-line retreats away from the obstructing wall in all simulations where it is present (see 
also lOka et al.ll2008l : lMurphyll2010l ). while generally staying near the lower base of the current sheet. 
The exception is case F with a strong guide fleld which shows that the X-line becomes located near 
the top exit of the current sheet. The X-line drifts towards the upstream r e gion with the stronge r 
magnetic field in most simulations with asymmetric infiow (e.g.. IUgaill2000l : ICassak &: Shavll2007l ). 
During simulations with symmetric upstream magnetic fields but asymmetric upstream densities, 
the X-line does not drift significantly along the inflow direction. 

For each simulation we find that the plasma fiow velocity at the X-line is substantially different 
from the time derivative of the X-l ine's position . Any differenc e between these two velocities must 
be due to resistive diffusion (e.g., ISeatonI l2008l : iMurphyl 120101 ). With the exception of the guide 
field and symmetric outfiow cases, the X-line retreats along the outfiow direction against the flow 
of the plasma. This occurs because diffusion of the infl ow component of the magnetic field along 
the infiow direction is able to shift the X-line position (JMurphvll20ld ). During reconnection with 
asymmetric inflow, the inflow component of plasma velocity at the X-line is in the same direction 
as the time derivative of the X-line position along the inflow direction. However, the plasma flow 
velocity is much faster than the rate of change in position of the X-line along the inflow direction. 
Consequently, one must be careful when converting simulations to the reference frame of the X-line 
because the plasma flow velocity will likely be different. 

While our simulations provide significant insight into the basic physics and observational con- 
sequences of line-tied asymmetric reconnection in the solar atmosphere, it is important to note 
the limitations of our models. Most noticeably, /3 is larger in the simulations than in most of the 
corona. Because the simulations start from a perturbed initial equilibrium, total pressure balance 
is enforced along the inflow direction by increasing plasma pressure on the weak magnetic fleld 
side. It is difficult to maintain a large asymmetry in low-/3 plasmas. The asymmetries in the simu- 
lations are probably more extreme than in the solar atmosphere. We speculate that magnetic fleld 
asymmetries of ~ 10-25% are common, with differences of perhaps ^50% possible in more extreme 
cases. The assumption of an initial equilibrium is also questionable since CMEs are non-equilibrium 
events. The simulations do not include radiative cooling, coronal heating (except for viscous and 
Ohmic heating), or initial vertical stratiflcation of the solar atmosphere. The outer conducting 
wall boundary conditions will affect the results on long time scales by additional line-tying as well 
as pileup of reconnection exhaust. The Lundquist numbers in our simulations are ^10^ so our 
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current sheets are below the S ^ 5 x 10 threshold for the onset of the plasmoid instability (e.g. 



Loureiro et al. 



2010 



2007 



Barta et al. 



20081 . 



2011 



Bhattachariee et al.1 '2009: Huang fc Bhattacharieell2010l: IShepherd fc Cassak 



Ni et all I201Q : iShen et al 



2011 



Biskamp 



198' 



ar 



Our simulations 



consequently have only one X-line. 



Future work on this problem should be performed using a combination of observations and 
improved numerical simulations. In particular, more realistic initial configurations will provide 
more detailed predictions that can be compared directly to observations. The HyLoop suite of codes 
(jWinterl l2009l : IWinter et al.ll201ll ) has the ability to take asymmetric loop configurations directly 
from these simulations, inject a population of energetic particles, and predict flare emissions in 
detail. Observationally, investigations of candle flame post-flare loops, current sheet drifting, and 
footpoint motion and relative intensities will provide constraints on the theoretical models. In 
particular, events that display multiple signatures will provide the most useful constraints and 
provide the most complete story. 

Many open questions remain for asymmetric reconnection in the solar atmosphere. They 
include: (1) How asymmetric are typical flare/CME current sheets? (2) What mechanisms are 
responsible for asymmetric HXR footpoint emission during flares? (3) What causes the drifting 
motion observed in many flare/CME current sheets? (4) Where is the principal X-line in these 
current sheets? (5) How important are CME current sheets to the eruption as a whole? (6) Are 
CME current sheets able to instigate circulation in the rising flux rope? (7) What are the effects of 
the patchy distribution of magnetic flux in the photosphere? We hope to address these problems 
in future work. 
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Table 1. Simulation Parameters 



Case 



VAh 



Ro /o -Bzo 



Notes 
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24 


80 


1.00 


1.0 
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60 


80 


0.71 


0.5 
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60 


80 


0.50 


0.25 
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60 


80 


0.35 


0.125 
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60 


80 


0.50 


0.25 
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0.25 
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32 
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0.63 


1.0 


4 





5040 
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60 


48 


0.50 


0.25 


1 





8000 



Symmetric inflow 

Asymmetry in the upstream magnetic field 
Asymmetry in the upstream magnetic field 
Asymmetry in the upstream magnetic field 
Like case C but with triple resistivity 
Like case C but with a guide field 
Symmetric magnetic field, asymmetric density 
Symmetric outflow, periodic in y-dircction 



^Using Equation I I15II assuming L - 
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Fig. 1. — Simulation results for case C with Rq = Blq/Brq = 0.25 at t = 100. Shown are (a) the 
magnetic flux, Az, (b) the out-of-plane current density, Jz, (c) the inflow component of velocity, 
Vx, (d) the outflow component of velocity, Vy, and (e) the plasma pressure, p. The 'o' denotes 
the position of the 0-point, and the 'x' marks the spot of the X-line. Only a portion of the 
computational domain is shown. 
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Fig. 2. — The position of the X-hne as a function of time for cases A-H along the (a) inflow and (b) 
outflow directions. Case H is not shown in the bottom panel because yn(t) = due to symmetry. 
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Fig. 3. — A comparison of the plasma flow velocity at the X-line, {Vx{xn,yn),yy{xn,yn)), to the 
time derivative of the position of the X-line, {dxn/dt,dyn/dt), for the (a) inflow direction and 
(b) outflow direction, for case C with Rq = 0.25. The plasma flow velocity at the X-line differs 
significantly from the rate of X-line motion along both directions. 
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Fig. 4. — Simulation parameters for a slice along the inflow direction at the position of the X-line 
for case C with i?o = 0.25 at t = 100. The vertical dotted line represents the position of the 
X-line along the inflow direction. Shown are (a) components of the out-of-plane electric field, (b) 
contributions to total pressure balance along the inflow direction, (c) components of the in-plane 
magnetic field, and (d) the inflow and outflow components of velocity. 
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Fig. 5. — Reconnection electric field strengths at the X-line as a function of time for cases A~H. 
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Fig. 6. — Shown are the magnetic flux (solid blue contours) and flow velocity vectors in the 
frame of the 0-point (red arrows) for case C at t = 100. The 0-point is denoted by the circle at 
{xo,yo) = (—1.21,14.53). The flow pattern about the 0-point is counterclockwise; consequently 
there is net vorticity in the outflow plasmoid. 
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(a) Case A, t=49.0 (b) Case B, t=45.5 (c) Case C, t = 58.1 (d) Case D, t = 99.4 




Fig. 7. — Magnetic flux contours for cases A-D when the X-hne in each simulation (denoted by 'x') 
is located at y = 2.9. The dashed green line gives the position of the loop apexes (where By = 0). 
The red dotted line represents the asymptotic prediction made by Equation (j25p . 
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Fig. 8. — Magnetic flux contours for the analytical potential field solution of asymmetric post-flare 
loops as given by Equation ()2ip very near the magnetic field reversal along the lower boundary 
at X = a (represented by the red dotted vertical line) . The dashed red line is the location of the 
loop-tops as given by Equation ([^S]) . The solutions are for magnetic field ratios of (a) i? = 1, (b) 
R = 0.5, (c) R = 0.25, and (d) R = 0.125. 



